APS/123-QED 



Unifying all mass discrepancies with one effective gravity law? 

HongSheng Zhao 

SUPA, School of Physics and Astronomy, University of St Andrews, KY16 9SS, Fife, UK 
Vrije Universiteit, De Boelelaan 1081, 1081 HV Amsterdam, Netherlands 



(N 

o 



(N 

o 

u 

6 



> 

(N 

m 
o 

(N 



Benoit Famaey 

Observatoire Astronomique, Universite de Strasbourg, CNRS UMR 7550, Strasbourg, France 

(Dated: July 27, 2012) 

A remarkably tight relation is observed between the Newtonian gravity sourced by the baryons and 
the actual gravity in galaxies of all sizes. This can be interpreted as the effect of a single effective force 
law depending on acceleration. This is however not the case in larger systems with much deeper 
potential wells, such as galaxy clusters. Here we explore the possibility of an effective force law 
reproducing mass discrepancies in all extragalactic systems when depending on both acceleration 
and the depth of the potential well. We exhibit, at least at a phenomenological level, one such 
possible construction in the classical gravitational potential theory. Interestingly, the framework, 
dubbed EMOND, is able to reproduce the observed mass discrepancies in both galaxies and galaxy 
clusters, and to produce multi-center systems with offsets between the peaks of gravity and the 
peaks of the baryonic distribution. 

PACS numbers: 98.10.+Z, 98.62.Dm, 95.35.+d, 95.30.Sf 



I. INTRODUCTION 

The dynamics of spiral galaxies exhibit a well- 
documented fine-tuned conspiracy between the surface 
density of baryonic matter and the total gravitational 
field |V$| [1]. Whatever the explanation for it, this obser- 
vational fine-tuned relation involves an acceleration scale 
ao ~ cHo/2ir ~ 10~ 10 ms~ 2 such that the effects of the 
putative dark matter appear at |V$| < ao and disappear 
at |V$| > ao in galaxies of all sizes. This is encapsulated 
by the empirical formula of Milgrom Q, at the basis of 
the Modified Newtonian Dynamics (MOND) paradigm. 

Nevertheless, in the central parts of galaxy clusters, 
where the observed acceleration |V$| > ao, the above 
prescription underpredicts the mass discrepancy by a fac- 
tor of a few. Putting aside the obvious possibility of 
this constituting a practical falsification of MOND, three 
other possibilities are that (i) there are dark baryons in 
the central parts of clusters, (ii) there is additional non- 
baryonic dark matter in the central parts of clusters, or 
(iii) Milgrom's formula is the limiting case of a more gen- 
eral empirical relation. For extensive discussions of (i) 
and (ii) see, for instance Here we rather concen- 

trate on the possibility of case (iii). 

A possibility to explain the discrepancy in galaxy clus- 
ters is that the effective force law needs a new scale in 
addition to ao. An obvious scale distinguishing galaxy 
clusters from galaxies is the deepness of the potential 
well |$|. For instance, Bekenstein Q recently proposed 
in this vein to add a velocity scale to Milgrom's prescrip- 
tion, such that the acceleration scale ao is a function of 
|$|. In this note, we explore this very general idea and 
digress on its observational consequences. We start by 
reviewing the MOND formalism in Sect. 2, then expose 
its generalization in Sect. 3. Plausible variations of ao as 
a function of $ and their observational consequences are 



then explored in Sect. 4 and 5 respectively, and conclu- 
sions are drawn in Sect. 6. 



II. MOND 

The idea of MOND is to link the observed gravita- 
tional attraction, g = — V$, to the Newtonian gravita- 
tional field calculated from the observed distribution of 
visible matter, g~N = — V</>jv, by means of an interpolat- 
ing function fj,: 
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where /jl(x) — > x for x <C 1 and fi(x) — > 1 for x ^> 1. 
However, this expression is not robust since it does not re- 
spect usual conservation laws. A consistent modification 
of gravity at the classical level should come from modi- 
fying the gravitational part of the Newtonian lagrangian 
density /^Newton = —p4>N — |V0at| 2 /(87tG). Bekenstein 
& Milgrom [lfj have developed a modified gravity frame- 
workwhere |V0at| 2 is replaced by a 2 l -F(|V$| 2 /a ! ), so that 
the modified lagrangian density reads: 
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where x — , and F can a priori be any dimensionless 
function. Varying the action with respect to $ then leads 
to a non-linear generalization of the Newtonian Poisson 
equation: 



p = V. 
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where [i(x) = F'{x 2 ). In order to recover the right lim- 
its for /J,(x), one needs to choose F{y) — > y for y 
1 and F(y) -> |y 3/2 for y < 1. 
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In the most popular covariant version of the theory 
[ilTfilj ] the weak-field potential is $ = <j) N + cj> where 4> 
is a scalar field whose lagrangian density is proportional 
to OqF (a^" 2 /i a/3 V Q 0V^0) , where h a @ is a combination 
of the metric and of a time-like unit vector. This scalar 
field lagrangian density is thus fully similar to Eq. [2j and 
any modification of the classical MOND lagrangian could 
thus immediately be translated into a modification of the 
relativistic version of the theory. 
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F(x 2 ), A = A (d>) 2 , (4) 



where x 2 
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and Aq(<&) plays the role of ao 



in the usual MOND. This leads to the following modi- 
fied Poisson equation, generalizing the MOND equation 
(Eq.©: 



p = V. 
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where fi(x) = F'(x 2 ), and Fi(y) = yF'(y) — F(y) — 
ydfi, which is typically a positive quantity of order 
unity. 

Interestingly, with this "extended MOND" (EMOND) 
formalism, Gauss' theorem (or Newton's second theorem) 
would no longer be valid in spherical symmetry. This 
means that masses outside a sphere can affect the local 



gravitational field strength, and the empirical EqQ] no 
longer holds rigorously even in spheres, due to the addi- 
tional term T2, although this term is often small as we 
shall see in the next sections. So we still have approxi- 
mately the EMOND counterpart of Milgrom's empirical 
law: 
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III. A MORE GENERAL FORCE LAW? 

Eqs. [5] and [3] suggest that the gravitational constant G 
is effectively a running function of the gradient modulus 
|V$| only. In other words, the lagrangian of the associ- 
ated scalar field </> has no potential and depends only on 
its kinetic term. Any direct dependency on the field $ 
would however be natural, since it is more general and, 
as we shall see, empirically desirable. 

As noted by Bekenstein [U , an obvious distinction be- 
tween galaxies and galaxy clusters is the deepness of 
the potential well. So one could imagine that the ef- 
fective transition-acceleration ao is not a constant but 
is a monotonically increasing function of the potential 
depth, e.g., exp ^/(lOOOkms" 1 ) 2 ! x 10 _10 ms~ 2 , which 
boosts the MOND effect in galaxy clusters by a fac- 
tor of a few[9|. However, a side effect of Bekenstein's 
exponentially-varying function is that it predicts a value 
of ao of the order of 10 10000 ms~ 2 once the potential iv. 
reaches the order of c 2 , i.e., a neutron star or a stellar 
black hole would exhibit an undesirable deep-MOND be- 
havior (see Fig.Q]). Its effect on large scale structures can 
be severe, too. It thus seems necessary to study a wider 
set of models for the effective variation of ao, and select 
the best one from a phenomcnological point of view. 

To extend the MOND lagrangian in this vein a very 
general way is to write the following lagrangian: 



Masses outside a spherical shell can also affect the lo- 
cal gravitational field strength by changing the bound- 
ary value of $ or the zero-point of the potential 
$00, hence changing Aq at the shell. The appear- 



ance of such an absolute zero-point is a curable ar- 
tifact of our introduction of a non-covariant Aq(&) 2 
instead of a covariant Aq((/)) 2 . Our lagrangian in 
Eq. 2] can be easily redesigned for a covariant scalar 
field, e.g., by replacing -A ($) 2 F(| V$| 2 /A ($) 2 ) 
-Al{4>)F (A {<j))~ 2 h afi V Q 0V (f). The boundary values 
of 4> and A ((j>) 2 are then well-defined by the cosmological 
evolution of (j)^. The zero point of the Newtonian po- 
tential (f>M can raise or lower 4>, but this has no effect on 
the scale A^ ((f)) or the acceleration or the light deflection 
power V<i> = V0 + V<pN- Note that we hereafter assume 
that the covariant version deflects light in a similar fash- 
ion as General Relativity, by invoking the same vector 
field as in TeVeS UM in addition to the scalar field 6 1221 



A REALISTIC EFFECTIVE VARIATION OF 
THE ACCELERATION SCALE 



Keeping the same philosophy as the a priori unknown 
/i-function of MOND in mind, we examine plausible po- 
tential dependencies of Aq(<&) from a phenomenological 
point of view. The only phenomenological requirements 
would be that (i) ^4o does not have a too large value 
in the strong field or cosmology regime, such that Gen- 
eral Relativity is effectively recovered there, (ii) Ao takes 
roughly the usual value of ao in galaxies and does not 
vary too much from galaxy to galaxy or inside a given 
galaxy, (iii) Aq is boosted in galaxy clusters by a rea- 
sonable amount compared to the usual value of ao so to 
make a real observational difference with MOND. 

Note also that isolated and non-isolated systems should 
be coherently transitioned from each other. This means 
that the potential $ (or scalar field 4>) in Ao should re- 
ally be the total potential due to all the material. This 
reminds of the external field effect (EFE) of MOND fl, 
but is more contriving since the external potentials can 
be non-neglibible even when the external force is negli- 
gible. So, the fact that one can neglect the weak EFE 
acting on some objects in MOND does not mean that 
one can ignore the external potential effect (EPE) of 
EMOND. The cosmological value of the scalar field in 
the covariant version would thus set a cut-off to Ao((j>) 2 . 
A crude estimate of the EPE is to consider that most 
galaxies and dwarfs are actually within some ~ 15 Mpc 
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FIG. 1: Top: shows the A (<E>) in units of 10~ 10 ms ~ 2 for sys- 
tems of increasing escape speed V cs c = y / 2 < f> — 2'i>oo for several 
possible functions: A (3>) = cH ^/ (2500) 2 ) (Eq.EJ for 
fi(x) = x(l + x n y 1/n with n = 2 (thick solid), n = 1,4, ...,oo 
(thin solid curves). Also shown are the truncated curve 
Ao($) = exp([|3>/(1000) 2 |, 2] min ) x 10" 10 ms- 2 [dashed] and 
the corresponding exponential untruncated curve [dotted]. 
Bottom: The value of fi(\\/^\/A (^)) vs. |V$| using Eq. [7] 
for systems with circular velocities of 1000 km/s (galaxy clus- 
ters, thickest curve) as well as 300, 200, 100 and 50 km/s 
(thin curves); the latter four are in the grey zone between 
the traditional forms (dotted curves) n{x) = x/(l + x) and 
fj,(x) = x/(l + x 2 ) 1 / 2 , acceptable in all galaxies. 



of a galaxy cluster (e.g., Virgo is a rich cluster ~ 15 Mpc 
from us). Such clusters will typically contribute some 
(300 km s -1 ) 2 to the potential at the boundary of galax- 
ies. So we can consider that the total |4>| is typically 
always above \<boo\ ~ 10 _6 c 2 , which will set Aq close 
to a lower cut-off with little variation for small internal 
potentials of small objects. Similarly, Aq is close to an 
upper cutoff for galaxies residing inside the deep poten- 
tial of a rich cluster, consistent with the deep-MOND like 
behavior of these cluster galaxies Q. 

While there are many plausible functional forms of 
Aq(&) 2 (see Fig. [TJ, a simple transition can be achieved 
by a /i-like dependency on the potential without intro- 



ducing any extra function: 

A = A ($) 2 = (cff )V(p) , p = |$c- 2 |/r 2 (7) 

with f3 ~ 8 x 10~ 3 . Here // is the same function as 
whichever /x-function of MOND, except that the argu- 
ment is the potential instead of the acceleration. To be 
specific we shall use fi(x) = x] (l + x 2 ) 1 / 2 hereafter unless 
otherwise stated. This empirical interpolating function 
simultaneously determines Aq through Eq. [7J The co- 
variant version could take A ((j)) 2 ~ (cH ) 2 [i (|</>/(/3c) 2 1) . 

With this choice, all the above requirements on a pos- 
sible realistic dependency of Aq on the potential are 
met. (i) Aq varies very little for large scale structure, 
and is almost equal to the observed amplitude of the 
cosmological constant. The value of Aq is limited to 
cHq ~ 8 x 10~ 10 ms~ 2 even for black holes, (ii) Due to 
the EPE described above, th e value of Aq does not go to 
zero but levels off at (cH ) VI$oo|/(/3c) 2 - lO^ms" 2 , 
i.e., the classical value of qq (see Fig.Q]) in small galaxies 
and in the solar system [14J. Additionally, there will be 
only very mild (< 50%) variations of Aq inside galaxies, 
as well as between galaxies, small variations which might 
be desirable to explain some anomalies in strong lensing 
galaxies fl5T4l9l ]. In any case this scatter is acceptable 
as the MOND formula works for all galaxies as long as 
an is within the range between 0.9 and 1.7 x 10~ 10 ms~ 2 
[HI, As illustrated on Fig. 1 (bottom panel), the ef- 
fective /i-function stays within the range allowed by data 
on galaxy scales, (iii) Finally, the much smaller fi on 
scales of galaxy clusters (bottom panel of Fig. Q} means 
a boost of MOND effect in these (see also next section). 

As for the Ti term in Eq. [5] further examinations show 
that T 2 = j§^F 2 , F 2 = pF"{p 2 ) (x 2 F'(x 2 ) - F(x 2 )) . 
In the deep-MOND regime we found F 2 < 0.1 for typical 
forms of n{p) and fi(x). So this extra density T 2 has 
a maximum of about 0.1/3 -2 ~ 1000 times above the 
cosmic critical density. Such a low over-density is not 
very significant even at large radii in galaxies, and is mild 
in galaxy clusters. 



V. BOOST OF GRAVITY AND OFFSETS IN 
GALAXY CLUSTERS 

In order to illustrate how such an effective gravity law 
would affect the dynamics of galaxy clusters, we concen- 
trate hereafter on the above /u-like dependency (Eq. [7]) . 
There are two basic questions to be answered: (i) how 
much dark matter would be inferred from a Newton- 
Einstein framework for a given distribution of baryons? 
(ii) Can this effective "phantom" dark matter be generi- 
cally offset from the baryons as observed in weak lensing 
maps of the Large Scale Structure and in objects such as 
the bullet cluster Q? 

In order to answer both questions, we consider a 
double peaked oblate potential: $(i?, z) = $oo + 
YXLl !500 2 In ^r 2 + 200 2 /2500, where n = Vz 2 + R 2 
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FIG. 2: Top: Displays for the multi-center galaxy cluster 
the density ratio of EMOND baryonic mass over Newtonian 
dynamical mass p har y 0n / ' Pdyna-mical ~ P ~ |V<E>|/A ~ 1/3 
(red, where the gravity peaks), ~ 1/30 (blue, where gravity 
is nearly zero). Bottom: Shows along the long axis Z, the 
acceleration scale Ao(3>) (thin solid), pdynamicai = V 2 <l?/47rG 
(dotted), and pharyon = Ti — T2 (thick solid, see Eq. 5) in 
10 6 M©kpc~ 3 . The ratio T 2 /T 1 (long-dashed) is shown to be 
small. 



and r-i = \J{z — 700kpc) 2 + B? . We then determine the 
dynamical density distribution V 2 $/47rG in Newtonian 
interpretation, and we determine the corresponding dis- 
tribution of baryons in EMOND with Eq. 5. The result 
is plotted on Fig. [5J the distribution of baryons is typi- 
cally about 7 times less massive than the dynamical mass 
inferred from the same potential in Newtonian gravity, 
and interestingly, this distribution of baryons is 4-peaked, 
with the main baryonic peaks clearly offset from the cor- 
responding peaks in the potential, or Newtonian dynam- 
ical mass peaks (see Fig.[2|). It is thus generally possible 



to create large "phantom" dark matter peaks generically 
offset from the baryon distribution within the present 
gravitational framework. This could be relevant to the 
modelling of the dark core in the trainwreck cluster [2l[ 
or to the bullet cluster [6j, provided that the covariant 
version deflects light in a similar fashion as General Rela- 
tivity, by e.g. invoking the same vector field as in TeVeS. 
We also note that the main baryonic peak seems to be 
generically in the more outer parts of the system com- 
pared to the main gravity peak, which is somewhat more 
similar to the situation of the train wreck cluster than to 
the bullet cluster. 



VI. CONCLUSIONS 

The effective gravity law explored in this note is able 
to boost gravity in galaxy clusters, eliminating the need 
for additional hot dark matter in MOND^a problem that 
has been known for more than 10 years [J, Q. Generally 
speaking, the new framework actually needs two interpo- 
lating functions, and two fundamental constants defining 
their transition domains. The framework is in principle 
very flexible, and the details of the interpolating func- 
tions can be narrowed (or excluded) by observational con- 
straints. It will for instance be important to test which 
forms can reconcile a large discrepancy in the clusters 
themselves with a not-too-large discrepancy in the clus- 
ter galaxies. Here, we recycled the standard /z function 
to tailor the A = Aq function (Eq.[7]): the result is a tran- 
sition between two plateaus, Aq ~ (cHq) 2 for large scale 
structures and Aq ~ a 2 , for small galaxies. To make the 
scheme (Eq. 2]) covariant and valid for lensing is straight- 
forward (see end of Section 3). Future detailed dynam- 
ical and lensing studies of galaxies and galaxy clusters 
should allow one to test whether this "extended MOND" 
(EMOND) scheme can indeed provide an effective gravity 
law reproducing mass discrepancies at all scales. 
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